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Abstract. Based on the Riesz definition of tlie fractional derivative tiie fractional 
Schrodinger equation with an infinite well potential is investigated. First it is shown 
analytically, that the solutions of the free fractional Schrodinger equation are not 
eigenfunctions, but good approximations for large k and for a ~ 2. The first lowest 
eigenfunctions are then calculated numerically and an approximate analytic formula 
for the level spectrum is derived. 
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1. Introduction 

Wave equations play a significant role in the description of the dynamic development 
of particles and fields; e.g. the Maxwell-equations describe the behavior of the electro- 
magnetic field in terms of coupled partial differential equations. In quantum mechanics 
a particle may be described by the non-relativistic Schr5dinger wave equation, where 
the kinetic term is given by the Laplace-operator. 

Fractional calculus introduces the concept of non-locality to arbitrary hitherto local 
operators [H] , [15] , [10] . This is a new property, which only recently attracted attention 
on a broader basis. The interest in a non-local dynamic description of e.g. quantum 
systems has been steadily increasing, because it is expected, that quantum phenomena 
may be treated more elegantly from a generalized point of view. 

Within this context, it is helpful to investigate fundamental properties of a fractional 
wave equation and to study general features of its solution. 

As an example, we will present in the following the main results of a solution of the 
fractional Schrodinger equation with infinite potential well. 
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2. The problem - analytic part 

Let the one dimensional fractional stationary Schrodinger equation in scaled canon- 
ical form be defined as 

- A"/2^(x) = (2.1) 

with the fractional Laplace-operator A"/'^. The definition of a fractional order deriv- 
ative is not unique, several definitions e.g. the Riemann[16j, Caputo[3], Liouville pT] . 
Riesz p!7j, Feller p], WevlfTS] coexist and are equally well suited for an extension of the 
standard derivative. 

In order to preserve Hermiticity for the fractional extension of the Laplace-operator |13j. 
we will explicitly apply the Riesz fractional derivative 

A"/'. -o-ZM = r(i + rn- + - y W + /(^ - Q,; (2.2) 

< a < 2 

where the left superscript in ^-D° emphasizes the fact, that the integral domain is the 
full space M. 

Since the eigenfunctions of the Riesz derivative operator are given as 

^zZ?" cos(A;x) = -|A:|"cos(A;x) (2.3) 

^2L>" sin(A;x) = -| sin(A;x) (2.4) 

the eigenfunctions of the potential free {V{x) = 0) fractional Schrodinger equation 
using the Riesz fractional derivative follow as 

Rz^+(A:,x) = cos(A:x) (2.5) 
Rz^-{k,x) = sm{kx) (2.6) 

where the it-sign indicates the parity and the corresponding continuous energy spec- 
trum follows as: 

= |A;|" k€R (2.7) 

Since the integrals dxcos{kx) and dxsm(kx) respectively are divergent, the 
eigenfunctions are not normalizable on the full domain M. 

In the case a = 2, which corresponds to the classical quantum mechanics, of course 
we may apply the classical box-normalization, which means we make a statement on 
the behavior of the wave-function isolated in a box of size e.g. 2q of the form: 

cos{kx)dx < oo (2-8) 
sin{kx)dx < oo (2-9) 

<? 

According to (2.3) and (2.4) a special feature of the Riesz derivative is that the type 



of the eigenfunctions of the free fractional Schrodinger equation does not change for 
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arbitrary a.. Hence the question arises, whether box-normahzation is a legal procedure 
in the case of non-local differential equations. At least, the physical meaning of box- 
normalization may be different. 

To avoid these interpretative difficulties we will switch on a potential such 
that the eigenfunctions vanish at infinity. An ideal candidate for such a potential is 
the infinite potential well, centered at the origin with finite size which is explicitly 
given by: 



V{x) = { - ^ (2.10) 




A reasonable ansatz for the corresponding eigenfunction follows as: 

\0 \x\>q 

Where the superscript ^^{x) emphasizes the fact, that the wave- function per defini- 
tionem is now confined inside the infinite potential well and vanishes outside. The 
superscript °°'I'(a:) indicates, that this function may at first be defined on the whole 
domain, but is used only inside the bounded domain of the infinite potential well. 
The normalization condition for ^'•^(x) is now given by: 



/+0O 
da;^*(x) = / dx°°*(x) = 1 (2.12) 
-oo J —q 



As a consequence, we may interpret ^'^(x) physically as a normalizable wave- function 

and its absolute value '^^'(a;)'-'^(x)* as a probability measure. 
In the classical, local case (a = 2) we obtain immediately 



U,T,± 



COS k^x \x\ < q and A; = 1, 3, 5, ... 

sink^x |x| <5and A; = 2,4,6, ... (2.13) 

\x\> q 



and the continuous energy spectrum changes to a discrete one, since A; G N is an integer 
now. 

ET' = iY/k^ k = 1,2,3,... (2.14) 

In order to solve the fractional Schrodinger equation of the infinite potential well for 
arbitrary a, we have to apply the Riesz derivative operator to ^'^{x). 
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For the positive semi-axis x > we obtain for x < q: 



TT 

^^+1 

sin(7ra/2) 

TT 



.sinvra/z r 
r(l + a) — ^ — '-^ X <^ 

/■g-^ oo^(^_g)_2°°vl/(x) + °°^(x + 
io 



a4 + 



q—x 



I q+x 

= "i^D^r^ix)) 0<x<q 

In addition, for x > q due to the non-local character of the Riesz derivative operator 
we obtain a non-vanishing finite term 

which is negligible only in the special case of the infinite potential well discussed here. 

The corresponding equations for x < just interchange the roles of — ^) and 

°°^{x + ^). Since parity is conserved for the infinite potential well, we may restrict to 



^,D-mx)) = r(l + ^ / X > q (2.16) 



(2.15) without loss of generality. 



We can write (2.15) in short-hand notation: 

^^DT^{x)) = \i,^D^r^{x)) (2.17) 

This may be interpreted as a modification of the Riesz-operator, which now only covers 
the inside region of the potential well. Obviously both operators differ significantly. 

~ i?" / ^.D- (2.18) 

This is a general feature of all fractional derivative definitions, which span over the full 
M e.g. Liouville's, Weyl's and Feller's definition. 

Only for the Riemann r-D" and the Caputo cD" definition of a fractional derivative 
and only for the very special case of an infinite potential well centered at the origin the 
operator is not altered: 

= ^D" (2.19) 



°°L>" = ^D" (2.20) 



c 

because the corresponding derivative definitions only cover the inside region of the 
potential well < x < q. 
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As a direct consequence, eigenfunctions of the free fractional Schrodinger equation 
based on the Riemann and Caputo derivative definition automatically determine the 
eigenfunctions of the same Schrodinger equation with infinite potential well fulfilling 
the additional constraint, that the functions should vanish at the boundaries of the 
infinite well [9j. 

These eigenfunctions ^^^^'^(x, a) are known analytically, normalizable on M for a < 
2 and are given in terms of the Mittag-LefHer functions Ea[z) and EaR{z) as: 



^^+(x,a) = xt-ii?„,.(-:E") (2.21) 

^,^-{x,a) = x'^-'E^^^i-x'^) (2.22) 

^vi/+(a;,a) = Eo^i-x'^) (2.23) 

'^^~{x,a) = a;t£;^i+|(-a;") 0<a<2 (2.24) 

which reduce to the trigonometric functions for a = 2 and the corresponding eigenval- 
ues are determined from the zeros of these functions [9], |10j . 

In contrast to this simple classical behavior we now consider the case of the fractional 
Schrodinger equation for the infinite potential well based on the Riesz definition of a 
fractional derivative. The following questions arise: 

• are plain waves still a solution for the infinite potential well? 

• if not, are they at least a good approximation? 

• what do the exact solutions look like? 

To answer these questions, we rearrange terms in the integral-operators In of the mod- 
ified Riesz-operator ^z-D": 

(nz^^r ^)) ix) = (r(l + a) ^^"(^"/^) (J, + + j3)oo^^ (^) (2.25) 
where In are given explicitely: 



^{x + - 2^'(a;) + ^(a; - 



rq+x _ ^) 

/3^(-) = / —^di = -2^{x) —di (2.28) 

J n—x S J q—x S 



^2^(^) = / ..+1 (2-27) 

J q—x S 
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Figure 1. The lowest solutions of —^y^D" cos{k-K / 2x) = Egk{a,x) for 
g = 1 in the range (dashed line) 2 > a > 0.25 (dotted line) in Aa = 0.25 
steps. For the classical, local case a = 2 the solution is indeed an 
eigenfunction, but for decreasing a deviations from cos(A;7r/2x) become 
more and more pronounced. On the other hand, for large k the error 
becomes smaller and cos{kTT /2x) becomes a good first guess for the exact 
eigenfunction. 

Using the test functions ^^{x) = cos(A;7r/2x) and ^~{x) = sin(fc7r/2x) we obtain 
with the help of [T], [5] and a cup of tea: 



h 



cos(fc7r/2x) 
sin(A;7r/2x) 



2a 



a - 2' L 
8 sin^(/c7r- 



1 



1^2 



a 3 a , o o ,q — X 

2-^;-A^V(^)2) + 



Aq 



2' 2' 2 
cos{k'ir/2x) 
sin(A;7r/2x) 



cos(A;7r/2x) 
sin(fc7r/2x) 

cos(A;7r/2x) 
sin(A:7r/2x) 



fcvr ^ fRe(T(x)) 
^2q' \-Im{T{x)) 



-2{q-x) 



a 



cos(A;7r/2x) 
sin(A;-7r/2x) 



(2.29) 
(2.30) 
(2.31) 
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Figure 2. The lowest numerically determined eigenfunctions ^'^(a, x) 
for q = 1 with positive parity for (dashed line) 2 > a > 0.25 (dotted 
line) in Aa = 0.25 steps. Only in the limit ^^(ce — )• 2,x) = cos(A;7r/2x) 
holds. 

with ii<2(a; b, c; x) is the hypergeometric function, Re(T(3;)) is the real part and Im(T(x)) 
is the imaginary part of the complex function T(x) given by 

T(x) = {-i)-e'"'^^^(ri-a,-ik7r^)-ri-a,-ik7r^)) (2.32) 

where T{a, z) denotes the incomplete T-function on the complex plane. 

The kinetic part of the fractional Schrodinger equation may therefore be calculated 
analytically and may be written with the testfunctions °°'^'^{x) = cos(/c7r/2x) and 
^^-(a;) = sin(/c7r/2x): 

-^2^"cos(A;7r/2x) = Egk{x) fc = 1,3,5,... (2.33) 

-^zD"sin(A;7r/2x) = Egk{x) A; = 2,4,6,... (2.34) 

We define a pseudo-normalization condition 

fffc(x = 0) = 1 A: = 1,3, 5,... (2.35) 

gk{x = l/k) = 1 A: = 2,4,6,... (2.36) 

which allows to compare gk{x) with the trigonometric functions. 
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Figure 3. The lowest solutions of — ^2-^°^ sin(A;7r/2x) = Egk{a,x) for 
g = 1 in the range (dashed line) 2 > a > 0.25 (dotted line) in Aa = 0.25 
steps. For the classical, local case a = 2 the solution is indeed an 
eigenfunction, but for decreasing a deviations from sin(/c7r/22;) become 
more and more pronounced. On the other hand, for large k the error 
becomes smaller and sin(A;7r/2x) becomes a good first guess for the exact 
eigenfunction. 



In figures [T] and [3] the results are sketched. The bad news is, that neither gk{x) = 
cos(fc7r/2x) nor gk{x) = sin(A;7r/2x) holds. Consequently the solutions of the free frac- 
tional Schrodinger equation based on the Riesz derivative definition are no eigenfunc- 
tions of the same fractional Schrodinger equation with infinite potential well. This 
answers the first question. 

The good news is, that the deviations become more and more negligible for k ^ 1. 
This answers the second question: For large k, in the vicinity of a ~ 2 and for x ~ 
the trigonometric functions seem to be a good first guess. 



As a consequence of the pseudo- normalization condition (2.35), we can give an ana- 
lytic expression for the approximate energy spectrum. With 



RZ 



D'^cosikTr/2x)\^=o = E^{a) A; = 1,3,5,... (2.37) 
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Figure 4. The lowest numerically determined eigenfunctions ^'^ (a, x) 
with negative parity for g = 1 in the range (dashed line) 2 > a > 0.25 
(dotted line) in Aa = 0.25 steps. Only in the limit ^^{a — t- 2,x) = 
sin(/c7r/2a;) holds. 



we obtain (since I2 vanishes for x = 0): 

E^ia) = yTia)i^{^,F4l-'l-,l.2-"^^^^^^ (2.38) 

which according to the above mentioned criteria will turn out to be a good approxima- 
tion of the exact energy eigenvalues for large k in the vicinity of a ~ 2. 



3. The solution - numerical part 




Since we have shown by explicit analytic evaluation of the fractional derivative, that 
the trigonometric functions are no eigenfunctions of the infinite potential well, but only 
good approximations, we will calculate the exact solutions numerically. 
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Figure 5. A comparison of the numerically determined energy values 



(triangles) with approximate formulas £'^(a) from (2.38) (solid line) 
and E^'^'^{a) from (2.7) (dashed line) normalized to q = 1. 



For that purpose, we expand the exact solution ^^{x) in a Taylor series 



N 



lim > 

Af-s>oo ^-^ 



a2nX 



„2n 



n=0 



2n+l 



(3.1) 
(3.2) 



n=0 



and insert it into the fractional Schrodinger equation for the infinite well potential: 

-^^D"^^{x) = E,^^{x) (3.3) 



The integrals on the left may be evaluated analytically and lead to transcendental 
functions in x, which then are expanded up to order in a Taylor-series too. This 
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leads to 



N N 

^62„x2" = i?fc5^a2„x2" (3.4) 

n=0 n=0 
N N 

j;62„+ix2"+i = i?,J^a2„+ix2"+i (3.5) 

n=0 n=0 

A term by term comparison results in a system of non-linear equations of the type 
{Ekttn = bn{an)} on the set of variables {an,Ek}, which is solved numerically. For 
practical calculations we set N = 20, which yields an accuracy of the calculated energy 
levels of about 0.25% for the ground state. 

Results are presented in figures [2]and|4} For a = 2 we obtain the classical trigono- 
metric functions. For decreasing a the eigenfunctions show a increasing tendency to 
shift towards the walls. This is exactly the behavior, which is not modeled by the 
trigonometric test functions presented in the last section. In figure [5] we compare the 
determined energy levels with the presented energy formulas. Especially E^{a) from 



(2.38) is a useful approximation. 

4. Conclusion 

We have demonstrated, that the non-local character of the fractional operators used 
in the fractional Schrodinger equation indeed needs special attention. Concepts like box 
normalization, WKB-approximation or piecewise solution may work well in a classical 
local approach, but cause errors when applied to non-local problems. 

On the other hand, we have shown, that such local strategies may lead to useful 
approximations e.g. in low-level (a = 2 — e) fractional problems. 

This may be understood from 

lim ii,Z)" = - (4.6) 

g— s>oo 

as a consequence of the specific weight w{^) = ^^^rr in the integral definition of the 
Riesz fractional derivative definition ( |2.2[ ) with the property 

lim w{0 = (4.7) 

The infinite potential well serves as a helpful tool to demonstrate the consequences 
of different approaches. 
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